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ABSTRACT 



We present an investigation of the morphology and arrival times of integrated radio 
pulses from the binary millisecond pulsar PSR J1022+1001. This pulsar is renowned 
for its poor timing properties, which have been postulated to originate from variability 
in its average pulse profile. Although a sub-class of long-period pulsars are known to 
exhibit mode changes that give rise to very large deviations in their integrated pro- 
files, this was the first millisecond pulsar thought to have an unstable mean profile. 
As part of a precision timing program at the Parkes radio telescope we observed this 
pulsar between January 2003 and March 2004 using a coherent de-dispersion system 
(CPSR2). A study of morphological variability during our brightest observations sug- 
gests that the pulse profile varies by at most a few percent, similar to the uncertainty 
in our calibration. Unlike previous authors, we find that this pulsar times extremely 
well. In five minute integrations of 64 MHz bands we obtain a weighted RMS residual 
of just 2.27 /iS. The reduced of our best fit is 1.43, which suggests that this pulsar 
can be timed to high accuracy with standard cross-correlation techniques. Combining 
relativistic constraints with the pulsar mass function and consideration of the Chan- 
drasekhar mass limit on the white dwarf companion, we can constrain the inclination 
angle of the system to lie within the range 37° < i < 56°. For reasonable pulsar 
masses, this suggests that the white dwarf is at least 0.9 M©. We also find evidence 
for secular evolution of the projected semi-major axis. 

Key v^fords: Pulsars, individual: J1G22+1001 



1 INTRODUCTION 

Pulsar timing is a highly versatile experimental technique 
that is used to provide estimates of everything from the 
age and magnet ic field strengths of pulsars, allowing pop- 
ulation studies llLvne et alJ Il985l: iLorimer et alJ ^93) , to 
micro-arc-second positio ns and tests of relativistic gravity 
l|van Straten et ail 1200 iT) . Central to the technique is the 
maxim that puls ars have average pu lse profiles that are in- 
herently stable jHelfand et al.lll975l) . which can be cross- 
correlat ed with a template to pro vide accurate times of 
arrival ([Manchester fc TavloJ [T977^ . Once corrections are 
made to account for the Earth's position with respect to 
the solar system barycentre, a model can be fit to the data, 
transforming pulsars from astrophysical curiosities into pow- 
erful probes of binary evolution and gravitational theories. 
One spectacul a,r example is the recently-discovered double 
pulsar system iLvne et al.ll20oJ) . 



The discove ry of the millisecond pulsars PSR B1937-I-21 
and B1855+09 feacker et alJll982l: Iseeelstein et al]ll98d) 
and their subsequent timing iKaspi et alJ Il994ft revealed 
that millisecond pulsars were very rotationally stable, rival- 
ing the best atomic clocks. Their lower magnetic fields ap- 
peared to give rise to a fractional stability that far exceeded 
that of the normal pulsar population (Arzoumanian et ^ 
Il994h . This stability, combined with their small rotation 
periods, gives the millisecond pulsars sub-microsecond tim- 
ing potential provided sufficient signal-to-noise ratio (S/N) 
pulse profiles can be recorded. This provides interest ing 
limits on a wide range of astrophysical phenomena iStairsI 
I2OO4I) . 

An array of accurately timed MSPs can even be used 
to constrain cosmological models and search for lo ng-period 
(nHz) gravitational waves jFoster fc BackeJll990D . Such an 
array requires frequent observations of a selection of pul- 
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sars that are known to time to a precision of less than 1 
/iS, preferably in two or more radio frequency bands. Imple- 
mentation of such an array would provide an opportunity to 
directly detect gravitational radiation in a regime that com- 
plements the sensitivity range of ground-based detectors, to 
test the accuracy of solar system ephemerides and to con- 
struct an entirely extra-terrestrial time scale. The question 
of whether a suitable sample of millisecond pulsars posi- 
tioned throughout the entir e celestial sphere ca n be selected 
from the pulsar catalogue jHobbs et alJl2004ll is under ac- 
tive assessment, through a collaboration between the Aus- 
tralia Telescope National Facility and Swinburne University 
of Technology. 

The feasibility of a timing array depends in part upon 
the intrinsic rotational stabilit y of the MSPs an d their 
lack of pulse profile variation. ICoenard fc Backed ll2004l) 
have recently discovered a "glitch" in the millisecond pul- 
sar PSR B1821-24, similar to those discovered in younger 
pulsars llShema r fc Lvnelll99(j) . Perhaps more disturbingly, 
lKramer'e^L^^999^ report that the binary millisecond pul- 
sar PSR J1022-f 1001 exhibits pulse shape variations that 
ruin its timing. Studies conducte d bv iKramer et al. (1999); 
iRamachandran fc KrameJ (j200^ argue that the pulsar ex- 
hibits significant changes in its pulse morphology on ~5 
minute time scales and narrow bandwidths. They inter- 
preted these variations as the source of the pulsar's unusu- 
ally poor timing prop erties. By mod e ling t he profile with 
Gaussian components, iKramer et alJ ^1999^ improved the 
timing and argued that the trailing component of the pulse 
was more stable than the leading component. 

The precision with which astronomers can predict pulse 
arrival times has been steadily improving over the past few 
decades with the advent of new technology and methodol- 
ogy. The charged interstellar medium disperses radio pulses 
and broadens the pulse profile across finite receiver band- 
widths. Astronomers crave both the wide bandwidths that 
permit high S/N integrated profiles, and the resolution of 
sharp features that permit strong cross-correlation with a 
standard template in order to minimise timing errors. Pul- 
sar astronomy has always benefited from adopting new tech- 
nologies that give increased time resolution and help defeat 
the deleterious effect of the interstellar medium on the prop- 
agation of radio pulses. 

Incoherent devices such as analogue filter-banks are 
prone to systematic errors due to interstellar scintil- 
lation and imperfect frequency responses. Digital spec- 
trometers can overcome some of these inadequacies but 
^Hankins fc Rickett ( 1975:) described a computational tech- 
nique (coherent de-dispersion) that removes pulse dispersion 
"perfectly" almost three decades ago. Unfortunately, at the 
time, computational hardware was not capable of process- 
ing all but the most modest of bandwidths. The exponential 
growth of computational power has permitted the develop- 
ment of new pulsar instruments that are capable of coher- 
ently de-dispersing large bandwidths in near real time. Such 
an instrument is CPSR2, a baseband recorder that permits 
near-real time coherent de-dispersion of 2x64 MHz bands 
using a cluster of high-end workstations IjBailcs 2003) . These 
new instruments represent the best opportunity to study 
small variations of integrated pulse profiles because they 
deliver an immunity against systematic errors induced by 
narrow-band scintillation and dispersion measure smearing. 



They also provide full polarimetry and multi-bit digitisation, 
radio frequency interference rejection and the opportunity to 
apply accurate statistical corrections that help elim inate sys- 
tematic errors induced by discretisation jjenet fc Anderson| 
Il998h . 

Motivated by the arrival of CPSR2, we have commenced 
a timing campaign of the best millisecond pulsars in an effort 
to make progress towards the aims of the timing array. PSR 
J1022-f 1001 was included because of its high flux density, 
and the paucity of better candidates at the hour angle it is 
visible from the Parkes radio telescope. PSR J1022-I-1001 is 
a recycled or millisecond pulsar with a pulse period of ap- 
proximately 16 ms. It orbits once every 7.8 days in a binary 
system with a companion whose minimum mass can be de- 
rived from the pulsar mass function, Eqn Q if we assume 
that the orbit is face on and the pulsar has a mass of 1.35 
M0. 

»m sin^ i 4it'^ sin^ i 



G 



(1) 



(mp -I- ttiwd)^ 

Here rrip and mwD refer to the mass of the pulsar and white 
dwarf respectively, i is the inclination angle of the orbit, Pb is 
the orbital period and G is Newton's gravitational constant. 
For PSR J1022-I-1001, this corresponds to a mass of at least 
0.7 Mq. 

Measurements of the dispersion measure along the line 
of site to the pulsar, combined with the Tavlor fc Corde^ 
1119931) galactic electron density model, place PSR 
J1022-I-1001 at a distance of roughly 600 pc from the sun, 
making it a relatively nearby source. Accurate timing of 
this pulsar could therefore lead to a greater knowledge of 
its orientation, companion masses, distance (via parallax) 
and proper motion. Unfortunately, this pulsar is very near 
the ecliptic plane and hence pulsar timing is only good at 
accurately constraining its position in one dimension. It is 
therefore a good target for very long baseline interferometry, 
as this could accurately determine the position and proper 
motion in ecliptic latitude. 

In this paper we present an analysis of the first 15 
months of CPSR2 observations of PSR J1022-I-1001, demon- 
strating that this pulsar can be timed to high accuracy using 
standard techniques. In section |5| we describe the observ- 
ing system and methodology, followed by a description of 
our data reduction scheme, including polarimetric calibra- 
tion. Section 13 describes the results of a sea rch within our 
data s et for variations of the type reported bv lKramer et alJ 
il999tl . Pulse arrival times calculated from this data set are 
analysed in section]^ followed by a summary of the newly 
derived pulsar spin and binary system parameters, including 
those that may become significant in the future. In section 
l^we discuss the implications of our work for this pulsar and 
precision timing programs in general. 



2 OBSERVATIONS 
2.1 Instrumentation 

The second Caltech Parkes Swinburne Recorder (CPSR2) is 
a baseband recording and online coherent de-dispersion sys- 
tem that accepts 4x64 MHz intermediate frequency bands 
(IPs). Using the in-built supercomputer, data can be pro- 
cessed either in real time if the pulsar has a dispersion 
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measure (DM) less than approximately 50 pc cm~'^ at a 
wavelength of 20 cm, or recorded to disk and processed of- 
fline. This allows a large amount of flexibility in observ- 
ing methodology. CPSR2 was commissioned in August 2002 
and has been recording data on a regular basis since about 
November 2002. Observing sessions are primarily conducted 
using a standard conflguration consisting of two indepen- 
dent 64 MHz bands, each with two orthogonal linear polar- 
isations. The down-conversion chain is configured to make 
both bands contiguous, at centre frequencies of 1341.0 and 
1405.0 MHz. In addition, a small number of observations 
are made at the widely separated frequencies of 3000.0 MHz 
and 685.0 MHz. After down-conversion and filtering to cre- 
ate band-limited signals, each IF is fed into the CPSR2 Fast 
Flexible Digitiser board which performs real, 2-bit Nyquist 
sampling. The digitisation process generates a total of 128 
million bytes of data every second. These samples are fed via 
Direct Memory Access (DMA) cards to two high-speed com- 
puters that divide the data into discrete blocks which are dis- 
tributed via gigabit ethernet to a cluster of client machines 
for reduction using a software program called PSRDISP 
ivan Stratcn 2003). PSRDISP implements software-based 
coherent de-dispersion faankins fc Rickettlll97§h . assuming 
prior knowledge of the DM along the line of sight to the pul- 
sar and the validity of the tenuous plasma dispersion law. 

2.2 Summary of Observations 

Observations were conducted at the Parkes radio telescope 
over a period of around 400 days between January 2003 and 
February 2004. The central beam of the Parkes Multibeam 
receiver dStavelev-Smith et alj[l993) was originally used at 
the front end, providing a system temperature of approxi- 
mately 21K at 20 cm. The Multibeam receiver was removed 
for maintenance in October 2003 and all 20 cm data since 
this date were recorded using the refurbished, wide-band H- 
OH receiver whose observable band encompasses that of the 
Multibeam receiver. The H-OH system is slightly less sen- 
sitive than the Multibeam, with a system temperature of 
25K at 20 cm wavelengths. In addition, a new coaxial dual- 
band 10/50 cm receiver system was installed in place of the 
Multibeam. During the 2004 observing sessions, data were 
taken with this system in order to expand our frequency 
coverage for the purposes of dispersion measure monitoring. 
Preliminary system temperature measurements of the coax- 
ial system yield values of 30 K and 40 K for the 10 cm and 
50 cm feeds respectively. 

Scheduled observing sessions were typically a few days 
in duration, occurring once or twice a month. Individual 
tracks of PSR J1022-I-1001 vary from approximately 15 min- 
utes to a few hours in duration. The data are somewhat bi- 
ased towards episodes of favourable scintillation as this al- 
lows the most efficient use of limited telescope time. In addi- 
tion to the on-source tracks, most observations were immedi- 
ately preceded or followed by a short (2 minute) observation, 
taken 1 degree south of the pulsar position, during which the 
receiver noise source was driven with a square wave at a fre- 
quency of 11.122976 Hz. These calibration tracks were used 
to characterise the polarimetric response of the signal chain 
so that corrections to the observed Stokes parameters could 
be made during data reduction. Observations of the radio 
galaxy 3C218 (Hydra A) were also taken at approximately 



monthly intervals (though often only at 20 cm) and used to 
calibrate the absolute flux scale of the observing system. 

2.3 Data Reduction 

Data from each dual-polarisation band were reduced to a co- 
herent filter bank with 128 spectral channels, each 0.5 MHz 
wide, corresponding to an effective sampling time of 2 /js. 
1024 phase bins were stored across the pulse profile, yield- 
ing a time resolution for PSR J1022-f 1001 of about 15 /^s 
per bin. Each individual data block recorded by CPSR2 rep- 
resents 16.7 seconds of data (1 Gigabyte of 2-bit samples). 
These blocks can later be combined to form integrated pro- 
files. For this purpose the auth ors u se the application set 
provided with the PSRCHIVE fHotan et al. 2004) scheme. 
Due to the lERS^ practice of retrospectively publishing cor- 
rections to the rotation rate of the Earth, the data presented 
in this paper are that for which corrections to the solar sys- 
tem barycentre could be accurately made at the time of 
preparation. Observations of PSR J1022-I-1001 are ongoing. 

2.4 Processing &: Calibration 

The set of 16.7 second PSR J1022-I-1001 integrations were 
summed to a total integration time of 5 minutes. These 
data were calibrated to account for the polarimetric re- 
sponse of the observing system, then all frequency chan- 
nels and Stokes parameters were combined to form a single 
total intensity profile corresponding to each 5 minute in- 
tegration. All of the receiver systems used during observa- 
tions of PSR J1022-I-1001 are equipped with orthogonal lin- 
ear feed probes, so the recorded polarisation products were 
corrected for relative gain and phase. Off-source calibrator 
observations were used to compute the relative gain and 
phase terms for each receiver system at the epoch of the 
observation, using a simple case of the scheme described by 
Ivan StraterJ ||200J). We do not attempt to correct for more 
subtle errors arising from cross-contamination between the 
two probes, primarily because sufficiently accurate models of 
all three receiver systems are not yet available. In the case of 
the central beam of the Multibeam receiver, existing models 
indicate that imperfections in the orthogonalit y of the two 
feed probes may be as large as several degrees ivan StraterJ 
2004). These imperfections break the fundamental assump- 
tion of orthogonality made when applying a simple complex 
gain correction and in duce err ors in the calculated total in- 
tensity of order 1-2 % JOrd et al. 2004|. Errors of this mag- 
nitude are only present in the proflle at phases where the 
fractional polarisation is close to unity. 



3 PULSE PROFILE STABILITY 

PSR J1022-I-1001 has an interesting profile morphology. At 
20 cm the profile consists of two sharp peaks (the princi- 
ple components), separated by about 0.05 phase units but 
joined by a bridge of emission (Fig . This characteristic 
do uble-peaked shape w as critical to the analysis performed 
bv lKramer et al.l (Il999ll . who make an effort to characterise 

^ http://www.iers.org 
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Figure 1. PSR J1022+1001 average profile at a wavelength of 20 
cm, formed from data taken with the central beam of the Parkes 
Multibeam receiver. The profile in the lower panel has been cal- 
ibrated for relative instrumental gain and phase. It has a centre 
frequency of 1341.0 MHz and covers a bandwidth of 64 MHz. The 
solid curve represents total intensity (Stokes I), the dashed curve 
represents linearly polarised emission and the dotted curve repre- 
sents circularly polarised emission. Note the high degree of linear 
polarisation in the trailing component. The relative (not abso- 
lute) position angle of the linearly polarised radiation is shown in 
the top panel. This profile has been calibrated against the flux of 
Hydra A (3C218). 

variations in the relative amplitude of the two components 
by normalising against one or the other and computing ra- 
tios of their height. The trailing pulse comp onent is observed 
to possess almost 100 % linear polarisat ion iOrd et al 
IStairs et aDllQOfll : Ixilotiris et alJll998h . whilst the leading 
component corresponds to a local minimum in the linear 
polarisation fraction which is near zero (Fig0. The pro- 
file also changes rapidly as a function of frequency, only the 
leading component is visible at a wavelength of 10 cm (Figl^J 
and the amplitude ratio at 50 cm (Fig|^ is nearly half that 
measured at 20 cm (Fig0. This spectral index variation 
within the profile and the highly asymmetric nature of the 
polarised fiux in the principle components, combined with 
the fact that PSR J1022+1001 is prone to strong episodes 
of scintillation, means that arrival time analysis over broad 
frequency bands is easily subject to the introduction of sys- 
tematic errors. It also complicates any search for variability 
by exacerbating the effect of instrumental imperfections, es- 
pecially those associated with polarimetry. 

3.1 Profile Normalisation 

In order to directly compare multiple integrated profiles 
from the same pulsar, brightness changes due to interstel- 
lar scintillation must be taken into account. This involves 
scaling or normalising each profile by finding a character- 
istic value associated with the source flux during the ob- 
servation and adjusting the measured amplitudes to ensure 
this value remains fixed across all observations. It should be 
noted that applying a simple relative scaling to all observed 
profiles masks any intrinsic variability in the total flux of the 
pulsar; but t hat is not the subj ect of this paper. Using the 
PSRCHIVE ( iHotan et alll2004|[) development environment. 
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Figure 2. PSR J1022+1001 average profile at a wavelength of 
10 cm, formed from data taken with the new coaxial 10/50 cm 
receiver system at Parkes. Again, the profile has been calibrated 
for relative gain and phase. It has a centre frequency of 3000.0 
MHz and covers a bandwidth of 64 MHz. Note the small degree of 
polarisation and the complete absence of the trailing pulse com- 
ponent. Flux calibrator observations were not made at the time 
these data were recorded, thus the amplitude scale is arbitrary. 
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Figure 3. PSR J1022+1001 average profile at a wavelength of 
50 cm, formed from data taken with the new coaxial 10/50 cm 
receiver system at Parkes. Again, the profile has been calibrated 
for relative gain and phase. It has a centre frequency of 685.0 
MHz and covers a bandwidth of 64 MHz. Note that the trailing 
component is almost 100 % linearly polarised and significantly 
stronger than the leading component. Flux calibrator observa- 
tions were not made at the time these data were recorded, thus 
the amplitude scale is arbitrary. 

software wa s constructed to perf orm tests similar to those 
described bv lKramer et alJ il999l) . In the process of analysis 
and testing, we developed a normalis ation scheme ba sed on 
the concept of diff e rence profiles (Hel fand et alJlT975^ . 

*Kra mer et alj 1^9^ describe a simple method of nor- 
malisation in which the characteristic quantity associated 
with each integration is taken to be the amplitude of one of 
the two component peaks. All the amplitudes in the profile 
are scaled by the constant factor required to give the cho- 
sen phase location the value of unity. In the case of PSR 
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J1022+1001, there is a choice as to which peak will be used 
as the reference. An alternative characteristic quantity is 
the total flux under the profile, or within a particular phase 
range. This quantity is simply the sum of all the profile am- 
phtudes within the region of interest. Normalising by fiux 
has the advantage that it uses no morphological informa- 
tion (other than the duty cycle of the pulsed emission in our 
case) and may therefore be better suited to detecting subtle 
profile variations. 

In this paper we compute fiuxes after subtracting a 
baseline level from each profile. A running mean with a 
phase width of 0.7 units is used to estimate the baseline 
flux, which is then subtracted from each bin. In order to 
further reduce our sensitivity to radio frequency interfer- 
ence and other factors that can distort the uniformity of 
an observed baseline, we normalise only by the flux in the 
on-pulse region of the profile. This region is defined by an 
edge detection algorithm that measures when the total flux 
crosses a threshold level. This edge detection is performed 
on the standard template profiles at each wavelength and 
the phase windows so defined are held fixed for the remain- 
der of the analysis. Unfortunately, normalising profiles by 
flux tends to artificially amplify noisy observations because 
a noise-dominated profile has a mean (and therefore a total 
flux after baseline removal) approaching zero. In the data 
reduction stage it is therefore advantageous to perform a 
cut on the basis of S/N. This also helps reduce contamina- 
tion by corrupted profiles. A S/N cut of 100 was deemed 
the best compromise between retaining a large fraction of 
the observed profiles and rejecting noise. This cut was ap- 
plied to the set of 5 minute CPSR2 integrations before any 
ana lysis was commenced . 

Ixramer et alJ il999l) 's technique of normalising to the 
trailing feature has a number of limitations: 

• It cannot place any limits on the variability of the fea- 
ture that is used to normalise the profiles. 

• Polarimetric calibration errors are maximised if the cho- 
sen component is highly polarised. 

• Finite signal to noise profiles will vary due to noise even 
if there is no intrinsic variability in the pulsar emission, this 
is true of all components in the pulse. 

We prefer our (more robust) flux-based normalisation 
scheme which normalises by the total flux in the on-pulse 
region. Flux normalisation allows us to properly identify 
which (if any) component is varying as it makes no prior 
assumptions. It also reduces the effect of poor calibration 
by averaging over the entire profile. There will still be a 
background of variability due to noise, but it will be spread 
uniformly. 

3.2 The Difference Profiles 

In order to detect subtle changes in pulse shape, it is possi- 
ble to subtract the amplitudes of a high S/N standard tem- 
plate profile from an appropriately scaled and aligned copy 
of a given integration. This procedure should yield Gaussian 
noise if the observed profile precisely matches the standard; 
any morphological deviations will protrude above the result- 
ing baseline. Before differences can be computed however, 
the observed profiles must be normalised to a standard tem- 
plate. We choose to normalise all profiles according to the 



r 0.0^5 




Pulse Phase 

Figure 4. The bottom panel shows superposed morphological 
difference profiles of 97 five minute total intensity (uncalibrated) 
integrations at a sky frequency of 1405.0 MHz. The middle panel 
shows the mean profile. The top panel shows the standard devia- 
tion of each phase bin (stars) in the difi'erence profiles displayed 
in the bottom panel and a series (filled circles) representing the 
mean value of the standard deviation in 64 windows, connected 
by lines, to aid perception of trends in the data. 



amount of flux in the on-pulse region, as described above. 
In addition, cross-correlation methods were used to deter- 
mine relative shifts, after which the observed proflles were 
rotated to align with the standard template. Because pro- 
file variability studies are sensitive to errors in polarimetric 
calibration, the difference profile analysis was performed on 
both the calibrated and uncalibrated profiles. Only the 20 
cm and 50 cm wavelengths had a sufficient number of high 
S/N profiles to make a difference profile study possible. 

Figs 0] - |7| show difference profiles computed from the 
CPSR2 observations at a wavelength of 20 cm, while Figs 
|S] & |U] show the 50 cm results. The number of calibrated 
profiles is less than or equal to the corresponding number 
of uncalibrated profiles, as not all observations have corre- 
sponding calibration tracks. 

As all figures show, the difference profiles are almost 
consistent with noise, with no alarming peaks near the lead- 
ing or trailing components. This can be seen from both the 
difference profiles and their standard deviations. The only 
evidence for variability comes from Fig |5] where ~2 % vari- 
ability above the noise floor can be seen near the trailing 
component. However, these variations are not present in 
Fig 2] which represents the same data set in the absence 
of polarimetric calibration. It would seem that the calibra- 
tion procedure itself is capable of inducing profile instability, 
which is perhaps not surprising considering the simplicity of 
the model used to correct for receiver imperfections. In ad- 
dition, CPSR2 has dynamic level setting that ensures the 
mean counts are equal in both polarisations. Provided the 
effective system temperature in the two polarisations is sim- 
ilar, calibration is almost unnecessary to form an accurate 
total intensity profile. It therefore seems a bizarre coinci- 
dence that these variations could be removed by simply ne- 
glecting to calibrate the data, and that they correspond to 
the maximum polarisation fraction in the profile. The sim- 
plest interpretation is that across 64 MHz bandwidths, this 
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Figure 5. As in Fig|3 superposed morpliological difference pro- 
files of 50 total intensity (calibrated) integrations at a sky fre- 
quency of 1405.0 MHz. The dashed line under the solid curve in 
the middle panel represents the total polarised emission in the 
mean profile. 
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Figure 6. As in Fig^ superposed morphological difference pro- 
files of 102 total intensity (uncalibrated) integrations at a sky 
frequency of 1341.0 MHz 



pulsar's mean profile is stable and that we are simply seeing 
the limitations of our calibration procedure. 



3.3 Peak Ratio Evolution 



iKramer et al.l 1^9^ infer the presence of smooth variations 
in the relative amplitudes of the principle components on 
time scales of a few minutes to an hour or more by calcu- 
lating peak amplitude ratios and demonstrating that they 
evolve at a level significantly above the uncertainty in the 
measurement. They also note that instances of such smooth 
variation may not be common and that the time scales in- 
volved can change from one data set to another. 

As a secon dary check, w e computed similar amplitude 
ratios to that of lKramer et a l. ( 1999) for our high S/N pro- 
files. A S/N cut of 30 was deemed sufficient for this analysis 
because it is less susceptible to baseline corruption than the 
difference profile test. The RMS of the off-pulse region in 



Figure 7. As in Fig|21 superposed morphological difference pro- 
files of 37 total intensity (calibrated) integrations at a sky fre- 
quency of 1341.0 MHz 




Figure 8. As in Fig|3 superposed morphological difference pro- 
files of 31 total intensity (uncalibrated) integrations at a sky fre- 
quency of 685.0 MHz 






Figure 9. As in Fig|5] superposed morphological difference pro- 
files of 31 total intensity (calibrated) integrations at a sky fre- 
quency of 685.0 MHz 
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Frequency (MHz) 


Ratio 


Std Dcv 


1405.0 


0.97 


0.07 


1341.0 


0.91 


0.07 


685.0 


0.58 


0.02 



Table 1. Ratio of leading to trailing component for the calibrated 
data set. 



Frequency (MHz) 


Ratio 


Std Dcv 


1405.0 


0.97 


0.06 


1341.0 


0.90 


0.06 


685.0 


0.58 


0.02 



Table 2. Ratio of leading to trailing component for the uncali- 
brated data set. 



904,02 904.04 904.05 904,08 904.1 904,12 

Time (MJD - 52000) 

Figure 10. Ratios of leading to trailing peak amplitudes during 
a single 2.5 hour observation taken on September 22, 2003, at a 
frequency of 1341.0 MHz. The error bars are computed based on 
the assumption that the amplitude of each peak could vary by 
the RMS of the noise in the baseline of the profile. 



each profile was used as an indication of the error in peak 
amplitude. The ratio of leading to trailing component ampli- 
tude was computed for all profiles and plotted against obser- 
vation time (see Fig llUl for an example). In addition, mean 
ratios and their associated standard deviation were com- 
puted over the 15 month time span of our data set. These 
values are presented in Tables 1 and 2. 

The evolution of the profile's components with fre- 
quency is quantified in Tables 1 and 2. It is interesting to 
note that calibration of the data increases the scatter in 
the measured ratios at higher frequencies, this further sup- 
po rts the notion that t he variability observed in Fig |S] and 
bv lKramer et al.1 il999l) is likely due to errors in the simple 
model used for polarimetric calibration. This does not seem 
to be the case at lower radio frequencies, as the calibration 
procedure has no detectable influence on the computed com- 
ponent ratios at 685.0 MHz. The 50 cm system may simply 
be better suited to the receiver model used for calibration, 
however the intrinsic scatter at this frequency is also much 
smaller so the effect of calibration may be imperceptible. To 
investigate the evolution of these component ratios in time, 
the data set was examined by eye in the hope of finding clear 
indications of non-random evolution between 5 minute in- 
tegrations. Fig 1101 shows the observation that exhibited the 
most variation. Calibration has little effect on these points. 
The variation about the mean is small and perhaps insignif- 
icant. 



3.4 Summary 

The CPSR2 data set suggests that across a bandwidth of 
64 MHz, PSR J1022-I-1001 does not vary significantly. Any 
instabilities, if present, must also be transitory in nature and 
therefore extremely difficult to characterise. It is possible 
that variability with a random or quasi-periodic structure, 
or small (< 64 MHz) characteristic bandwidth might still 
be present, but if so it would appear to have little effect 
on the profile when all frequency channels are combined. 
Profile instabilities should also be reflected in the timing of 
the pulsar, which we investigate in the next section. 



4 ARRIVAL TIME ANALYSIS 

iKramer et alJ (^9^ report that arrival times recorded at 
the Arecibo, Jodrell Bank and Effelsberg radio telescopes 
yield a best fit RMS residual of 15-20 us when applied to a 
model of PSR J1022-I-1001 and the binary system in which 
it resides. To achieve this RMS residual, a complicated stan- 
dard template profile consisting of five Gaussian components 
with fioating amplitudes was used, to compensate for sup- 
posed intrinsic profile variability. In contrast, the CPSR2 
data set yields arrival times (from each 5 minute integra- 
tion, across multiple radio frequencies) that fit our timing 
model with a RMS residual of 2.27 pis. These times of arrival 
(TOAs) were obtained using a Fourier domain cross correla- 
tion process. Simple, static standard template proflles were 
constructed from the sum total of multiple integrations and 
the baseline noise was flattened to zero, to decrease spu- 
rious self-correlation of the timing profiles with the noise 
in the standard, which can lead to an RMS that is artifi- 
cially low. Separate standards were used for each frequency 
band and were aligned to a common fiducial point. Our 5 
minute (uncalibrated) integrations were fit to a model of 
PSR J1022+1001 using the standard pulsar timing package 
TEMPO^. Fig llll shows the corresponding timing residuals. 
It is interesting to note that although the front-end receiver 
system was changed mid-way through the data set, there 
does not seem to be any large systematic offset between the 
two receivers and no jumps were used to fit across the bound- 
ary. Changes in cable length or amplifier response between 
the two systems must have some impact on the assignment 
of arrival times, however it would appear that the offset is 
too small to measure in this data set. 

A binary mode l of the type described by 
iBlandford fc Teukolskvl I^^T^ was used to model the 
spin-down characteristics of PSR J1022-I-1001 and the 
perturbations introduced by its companion. A list of fitted 
parameters and their corresponding values (with error 
estimates) is presented below (Table 3). The errors in each 
TO A were uniformly scaled by a factor of 1.195 before 

^ http:/ /www. atnf.csiro.au/research/pulsar/tempo/ 
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Figure 11. Timing residuals from all 5 minute integrations with 
S/N in excess of 30 and timing error less than 7 /xs. The RMS 
residual is 2.27 /xs and the reduced of the weighted fit is 1.43. 
The vortical scale i s an order of magnit ude better than arrival 
times presented by iKramer et alj il999t) who saw residuals as 
large as 400 /xs. 



fitting, to ensure the reduced per degree of freedom was 
equal to unity. This is necessary to produce better error 
estimates for the fitted parameters. Sky position is shown 
in both Ecliptic and Equatorial (J2000) coordinates. PSR 
J1022+1001 lies in the ecliptic plane, making it difficult 
to determine the ecliptic latitude (/3) accurately through 
timing measurements. This accounts for the relatively high 
error in our measurement of /3 in Table 3. In addition, the 
CPSR2 data has a time baseline of onl y 15 m onths, whereas 
the data presented by iKramer et al.l (|l99flh extends over 
more than 4 years. As a consequence, we do not fit for a 
proper motion in our model, choosing instead to hold this 
parameter fixed at t he value of -17 ± 2 mas yr~^ quoted by 
IKramer et all |l99i). With proper motion thus constrained, 
we obtain a significant estimate for the parallax of the 
system. Given the short temporal baseline of our data set, 
this parallax estimate should be considered preliminary and 
it is included because it appears significant. 

Fig im shows that PSR J1022+1001 has benefited from 
more regular observations in recent months. There are only 
a few small groups of points in the earliest part of the data 
set. This is another reason for not including proper motion 
in the model, as the small number of points observed during 
the beginning of 2003 would unreasonably dominate the fit. 
Parallax measurements are sensitive to day of year coverage 
more than total observing time, so it is still possible that the 
value of 3.3 ± 0.8 mas is meaningful. It is interesting to note 
that our parallax measurement implies a distance of only 
300 lfin° PC, as opp osed to the value of 600 pc derived from 
the iTavlor fc Cord es ( 1993) electron density model. Many 
pulsars within 1 kpc of the Sun are inconsistent with the 
fTavlor fc Corded il993i'l model at this level. If the system is 
indeed only 300 pc from the sun, other secular changes may 
soon be detectable in the orbital parameters, such as asini. 

The accuracy with which we have been able to time 
this pulsar affords the chance to derive new limits on sev- 



Parameter 


Value 


T^lplinf IP T ,rtT\ \ \ \ iHpo"! 

i_J 1.^11 LI LIL^ i^Wli* \ / \ ^^^^ / 




Rrlintir T,a1- ( fT\ Mpp-^ 


_n nfi3Qi ffi"! 


JT HJUt-l IVlULUJll 111 y\ t llldt) VI 1 


1 7 (9\ * 


Parallax (^rQas^ 


O.O \0 ) 


Period (ms) 


16.4529296931296 (5) 


JT Cl lUtl J-^cl 1 VdLl V t- v-LvV } 


4.33 (1) 


Period Epoch (MJD) 


52900 


. / _3 > 
Ulspersion IVieasure 1 cm pc 1 


-LU. ZuJ-OU 1 ) 


Projected Semi-Major Axis (It-s) 


16.7654148 (2) 


Eccentricity 


0.00009725 (3) 


Time of Periastron Passage (MJD) 




Orbital Period (days) 


I .OUOloUlOU 


Angle of Periastron (deg) 


97.73 (1) 


Right Ascension (a) 


10:22:58.015 (5) 


Declination {&) 


+10:01:53.2 (2) 


Number of TOAS 


555 


Total 


545.74 


RMS Timing Residual (/xs) 


2.27 


MJD of first TOA 


52649 


MJD of last TOA 


53109 


Total Time Span (days) 


460 



Table 3. PSR J1022+1001 iBlandford Teukolskvl il97fl) tim- 
ing model parameters derived from 15 months of CPSR2 observa- 
tions. The error in the last significant digit is given in parenthese s 
after the value. (* This value was given bv lKramer et alJ 0113)) 

eral physical parameters. Of particular interest is the possi- 
ble presence of orbitally modulated Shapiro delay as the dis- 
tance each pulse must travel into the companion star's grav- 
ity well changes throughout the 7.8 day orbit. Whilst fitting 
for the Shapiro delay parameters mwD and i is not directly 
possible due to the small amplitude of their timing signature, 
a statistical analysis of the effect these parameters have on 
the of the fit can still rev eal important information. The 
parameters derived from our IBlandford fc Teukolskvl ill 9 7d) 
model were incorporated into a model JPamou^^DerueUa 
^985, 1986) that includes the two Shapiro delay parameters. 
This new model was used to create a map for varying 
companion masses and inclination angles (Fig |12ll . This map 
allows us to place upper and lower bounds on the likely incli- 
nation angle, although it does not tightly constrain the com- 
panion mass. Based on 2-cr Shapiro delay contours, cosi > 
0.56. Assuming a neutron star mass of 1.35 M0 and a white 
dwarf companion (a valid assumption given the small eccen- 
tricity of the system), the Chandrasekhar mass limit for the 
companion constrains cosi < 0.8. These two limits combined 
imply that the inclination angle lies between 37° and 56° and 
suggest that mwD > 0.9 Mq. This suggests that the initial 
binary was near the limit required to produce two neutron 
stars, and that the white dwarf is composed of heavier ele- 
ments (maybe even ONeMg) than many of the millisecond 
pulsar companions which are most often He white dwarfs. 

T he model parameters presented in IKramer et alJ 
il999l) include a value for the projected semi-major axis of 
the orbit, x = asini = 16.765409(2) s. Given that the ref- 
erence epoch of our data set is some 7 years ahead of the 
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Figure 12. Shapiro delay map. Tlie thiclc lines represent 1-a 
contours, followed by 2,3 and 4-(t contours on tlie left hand side. 
The dashed line represents the mass function constraint, assuming 
a 1.35 M0 neutron star. 

corresponding iKramer et alj il999l) epoch, it is possible to 
compare our value of x in the hope of detecting a significant 
change as the proper motion of the system alters our line 
of sig ht to the plane of th e orbit. We find x = 0.8 ± 0.3 /^s 
vr' MSandhu et alJ (Il997l) show that: 

X = X coti{~ iiaSinQ, ~\- us cosQ,)ss~^ (2) 

where fia and represent the components of proper motion 
in RA and Dec and is the longitude of the ascending node. 
On ly the compon e nt of i i in ecliptic longitude was measured 
by iKramer et al.l lll999ll , but given our constraints on the 
inclination angle i, the only two unknowns in the equation 
are now and . If VLBI measurements could provide a 
value for /ia (and perhaps confirm our detection of parallax) , 
we would be able to place limits on the angle fi, further 
constraining the 3-dimensional orientation of the orbit on 
the sky. 

The reduced of 1.43 indicates that the arrival times 
are relatively free from unmodeled systematic effects. In- 
creasing the error estimates by 20 % gives a of unity, 
which is low by precision pulsar timing standards. Another 
way to view this low level of systematics is to plot the resid- 
ual against the arrival time measurement uncertainty for 
each TOA fFig |13ll . With perhaps one exception, there are 
alm ost no points many s tandard deviations from zero. 

Kramer et alJ lll999h argue that the leading component 
of PSR J1022-I-1001 varies across a characteristic bandwidth 
of approximately 8 MHz. In order to test this, a separate tim- 
ing analysis was performed on the 20 cm wavelength data 
before complete summation of the frequency channels, re- 
ducing the bandwidth in each integrated profile to 8 MHz. If 
the pulsar is varying on these bandwidths and 5 minute time 
scales, this should lead to an enormous increase in the RMS 
residual that grossly exceeds the expected factor of \/8. We 
find that the RMS residual of these 8 MHz bandwidth ar- 
rival times is 6.65 /is, which is very close to -\/8x 2.27 = 6.42 
/^s, implying that the increased scatter is consistent with 
the loss of S/N associated with the reduction in bandwidth. 
This analysis is strong evidence that there are no profile 
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Figure 13. Model residual plotted as a function of error in arrival 
time. 

variations with a characteristic bandwidth of approximately 
8 MHz and time scales of a few minutes. 



5 CONCLUSIONS 

PSR J1022-I-1001 has an interesting pulse profile morphol- 
ogy in terms of characteristic shape, polarimetric structure 
and evolution with radio frequency. These factors conspire 
to make it a difficult source to calibrate and analyse. Ex- 
tensive studies of the morphological differences between 5- 
minute integrations observed with CPSR2 at the Parkes 64 
m radio telescope are consistent with a stable pulse profile. 
This bodes well for the future of precision timing of both this 
source and millisecond pulsars in general. Although the RMS 
residual presented is probably not small enough to make an 
immediate contribution to any pulsar timing array, length- 
ier integrations and continued monitoring may yet push the 
timing of this pulsar below the 1 /is mark. 

We have used our improved timing to place some inter- 
esting limits on the geometry of this source which demon- 
strate that the inclination angle lies within the range 37° < 
i < 56°. If our parallax is confirmed, this system will make 
a good target for VLBI observations which could both im- 
prove the distance to the pulsar and determine the as yet 
unknown component of proper motion in the direction of 
ecliptic latitude. This would provide additional limits on the 
three dimensional orientation of the orbit via consideration 
of the change in projected semi-major axis. In the future, 
limits on x will improve and on a 10-year baseline, other 
relativistic observables may become measurable. We antici- 
pate that with another 20 months of timing, we should be 
able to obtain a more meaningful parallax and independent 
proper motion for this source in right ascension. Our error 
in uj is only 0.01°, which suggests that the expected rate of 
advance of periastron for this source u) ~ 0.01°yr~^ may be 
measurable on a 10-year baseline. 
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